Momentum spectroscopy of ID phase fluctuations in Bose-Einstein condensates 
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We measure the axial momentum distribution of Bose-Einstein condensates with an aspect ratio 
of 152 using Bragg spectroscopy. We observe the Lorentzian momentum distribution characteristic 
of one- dimensional phase fluctuations. The temperature dependence of the width of this distribution 
provides a quantitative test of quasi-condensate theory. In addition, we observe a condensate length 
consistent with the absence of density fluctuations, even when phase fluctuations are large. 
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One of the most striking features of Bose-Einstein con- 
densates is their phase coherence. Extensive experimen- 
tal work on dilute atomic gases has demonstrated a uni- 
form phase of three-dimensional (3D) trapped conden- 
sates p], 0, even at finite temperature |jg. In low di- 
mensional systems, however, phase fluctuations of the or- 
der parameter are expected to destroy off-diagonal long 
range order (see 0, 0] and references therein) . This phe- 
nomenon also occurs in sufficiently anisotropic 3D sam- 
ples, where phase coherence across the axis (long dimen- 
sion) is established only below a temperature that can 
be much lower than the critical temperature T c [(j . In the 
range < T < T c , the cloud is a "quasi-condensate", 
whose incomplete phase coherence is due to thermal exci- 
tations of ID axial modes, with wavelengths larger than 
its radial size. Quasi-condensates in elongated traps have 
been observed by Dettmer et al. |Jj, who measured the 
conversion, during free expansion, of the phase fluctu- 
ations into ripples in the density profile. Although the 
conversion dynamics is well understood || , the measured 
amplitude of density ripples was a factor of two smaller 
than expected. 

In this Letter, we report on the measurement of the 
axial coherence properties of quasi-condensates via mo- 
mentum Bragg spectroscopy. In previous work using 
this technique 0, 0, the finite size and mean- field en- 
ergy were the primary contributors to the spectral width. 
By contrast, the dominant broadening in our condi- 
tions results from thermally driven fluctuations of the 
phase, which reduce the coherence length 0, [(J. In- 
deed, the axial momentum distribution is the Fourier 
transform of the spatial correlation function C(s) — 
J d 3 r(*t( r -su z /2)$(r + su z /2)) |Tr|, where u z is the 
axial unit vector. When phase fluctuations dominate (i.e. 
T 3> T^), the axial momentum width is hence propor- 
tional to h/Lfi, where is the characteristic axial de- 
cay length of C(s). Experimentally, for 6 < T/T^ < 36, 
we find momentum distributions with Lorentzian shapes, 
whose widths increase with T. Such a shape is charac- 
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teristic of large phase fluctuations in ID which re- 
sult in a nearly exponential decay of C(s). Moreover, 
the momentum width agrees quantitatively with theo- 
retical predictions to within our 15% experimental un- 
certainty. This implies, in this temperature range, a 
coherence length substantially smaller than the quasi- 
condensate length 2L, from about L/18 to L/4. We have 
also checked an important feature of quasi-condensates: 
the suppression of density fluctuations even in the pres- 
ence of large phase fluctuations. 

We produce a Bose-Einstein condensate of 87 Rb in 
the 55*1/2 \F = L m F = —1) state. A new design of 
our iron-core electromagnet, with respect to our previ- 
ous work 0], allows us to lower the bias field of the 
Ioffe-Pritchard trap to obtain tighter radial confinement. 
Final radial and axial trap frequencies are respectively 
Lu ± /2ir = 760(20) Hz and u z /2ir = 5.00(5) Hz. The 
condensates, containing around 5 x 10 4 atoms 0], are 
needle-shaped, with a typical half-length L ~ 130 
and radius R ~ 0.8 /an. The chemical potential being a 
few times hw±, the clouds are between the 3D and ID 
Thomas- Fermi (TF) regime 0|. However, the low- lying 
excitations of the condensate are ID in character, due to 
the large aspect ratio of the trap [Til ]. 

Our momentum distribution measurement is based on 
four-photon velocity-selective Bragg diffraction. Atoms 
are diffracted out of the condensate by interaction 
with a moving standing wave, formed by two counter- 
propagating laser beams with a relative detuning 5 0, E| • 
Due to the Doppler effect, the momentum component res- 
onantly diffracted out of the condensate is p z = M(5 — 
8w R )/(2fc L ) with lu r = hkl/(2M), M the atomic mass, 
and /cl = 27r/A (A = 780.02 nm). The lasers are tuned 
6.6 GHz below resonance to avoid Rayleigh scattering. 
The laser intensities (about 2 mW/cm 2 ) are adjusted to 
keep the diffraction efficiency below 20%. 

To build the momentum spectrum of the quasi- 
condensate, we measure the fraction of diffracted atoms 
versus the detuning S between the counter-propagating 
laser beams. The differential frequency 5 must be sta- 
ble to better than the desired spectral resolution, about 
200 Hz for our typical = 10 ^m. The optical setup is 
as follows. A single laser beam is spatially filtered by a 
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fiber optic, separated into two arms with orthogonal po- 
larizations, frequency shifted by two independent 80 MHz 
acousto-optic modulators, and recombined. The modula- 
tors are driven by two synthesizers stable to better than 
1 Hz over the typical acquisition time of a spectrum. The 
overlapping, recombined beams are then sent through the 
vacuum cell, parallel (to within 1 mrad) to the long axis 
of the trap, and retro-reflected to obtain two standing 
waves with orthogonal polarizations, moving in opposite 
directions. After we mounted the critical retro-reflecting 
mirror on a long, rigid plate to minimize axial vibrations, 
the average over ten beat notes had a half-width at half- 
maximum (HWHM) of 216(10) Hz for a 2-ms pulse 0. 

The following experimental procedure is used to ac- 
quire a momentum spectrum. At the end of forced evap- 
orative cooling, the radio frequency knife is held fixed for 
6.5 s to allow the cloud to relax to equilibrium. Indeed, 
we observe axial shape oscillations of the cloud, trig- 
gered by the onset of Bose-Einstein condensation |17lll8| . 
despite a slow evaporation (less than lOOkHz/s) across 
T c . The magnetic trap is then switched off abruptly, in 
roughly 100 /is, and the cloud expands for 2 ms before 
the Bragg lasers are pulsed on for 2 ms. Since there are 
two laser standing waves moving in opposite directions, 
two diffracted components are extracted. We wait for a 
further 20 ms to let the diffracted atoms separate from 
the parent condensate, and take an absorption image 
(Fig. nji). Diffraction efficiency is defined as the frac- 
tion of atoms in each secondary cloud. We repeat this 
complete sequence for several detunings (typically 15), 
several times (typically 5). After averaging the diffrac- 
tion efficiencies measured at each detuning <5, we obtain 
two "elementary spectra", one for each diffraction com- 
ponent. 

We take the Bragg spectrum after expansion rather 
than in the trap to overcome two severe problems. In 
the trapped condensate, first, the mean free path (about 
10 ^m) is much smaller than its axial size, typically 
260 nm, so that fast Bragg-diffracted atoms would scat- 
ter against the cloud at rest |l9f. Second, the inhomo- 
geneous mean field broadening [2| would be of the order 
of 300 Hz, i.e. larger than the spectral width expected 
from phase fluctuations. By contrast, after 2 ms of free 
expansion, the peak density has dropped by two orders 
of magnitude |20|. and both effects become negligible. 
The measured momentum spectra faithfully reflect the 
phase fluctuations of the trapped quasi-condensate. In- 
deed, the phase fluctuations do not significantly evolve 
in 2 ms, since the typical timescale for their complete 
conversion into density ripples varies from 400 ms to 15 s 
for the range of temperatures we explore [111. Also, the 
mean field energy is released almost entirely in the radial 
direction, because of the large aspect ratio of the trap 
[itij l , and contributes only about 50 Hz of Doppler broad- 
ening in the axial direction. The only perturbation due 
to the trap release seems to be small axial velocity shifts 
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FIG. 1: (a) Absorption image of degenerate cloud (center) 
and diffracted atoms (left and right), averaged over several 
shots, after free-flight expansion, (b) Diffraction efficiency 
(rescaled) versus relative detuning of the Bragg lasers at 
T = 261(13) nK, corresponding to T/T$ = 20(2). A typical 
statistical error bar is shown. This spectrum is the superposi- 
tion of 12 "elementary spectra", as described in the text. The 
true average center is 30.18(2) kHz, close to 30.17 kHz, the 
four-photon recoil frequency. The solid line is a Lorentzian fit, 
giving a half-width at half-maximum (HWHM) of 316(10) Hz. 
(c) and (d) show respectively the residual of a Lorentzian and 
of a Gaussian fit to the above spectrum. Residuals are folded 
around the 8 = axis, and smoothed with a six-point-wide 
sliding average. 



(around 100/im/s) attributed to stray magnetic gradi- 
ents that merely displace the spectra centers. 

Bragg spectra have been taken at various temperatures 
between 90(10) nK and 350(20) nK, while T c varied from 
280(15) nK to 380(20) nK. The temperature was fixed to 
within 20 nK by controlling the final trap depth to a pre- 
cision of 2 kHz, and measured from a fit to the wings of 
an averaged absorption image. The fitting function is 
an ideal Bose distribution with zero chemical potential, 
plus an inverted parabolic profile for the quasi-condensed 
cloud. At each temperature, pairs of elementary spectra 
(described above) were collected across a 125-ms-wide 
range of hold times to average over residual oscillation 
and slowly varying fluctuations. All elementary spectra 
corresponding to the same temperature are reduced to 
the same surface, background, and center, and super- 
posed (Fig.Q]b). 

The line shape of the resulting spectra is clearly closer 
to a Lorentzian than to a Gaussian (see Fig. This 
is a significant result, because a Lorentzian-like profile 
is expected for a momentum distribution dominated by 
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phase fluctuations (see an d below), in contrast to 
the gaussian-like profile expected for a pure condensate 
From the Lorentzian fit, we extract the measured 

half-width Avyi for each temperature^ 

The theoretical results obtained in [(J yjj do not apply 
directly to our experiment. In those works, a 3D Thomas- 
Fermi density profile was assumed, whereas in our case 
modification of the density profile by the thermal cloud 
and radial quantum pressure must be accounted for. We 
find (and discuss below) that a parabolic profile is still 
a good fit function, but that the usual T = relations 
between n, L, R, and the number of condensed atoms Nq 
are no longer valid. We therefore extend the calculation 
of the axial correlation function in |ll| to an arbitrary 
density profile, in theJocal density approximation. In 
the mean field regime jj], from the result of [2l| for a ID 
uniform Bose gas at finite temperature, we obtain: 



C{s) 



dz ni(z) exp( 



2ni(z)L ct> / 



(1) 



where n\(z) = J d 2 r± no(rj_, z) is the axial ID density of 
the quasi-condensate, while no(r) is its 3D density pro- 
file. The coherence length near the center of the trap is 
given bYjL> = ft 2 ni(0)/(Mfc B T) y££. Following Petrov 
et al. U,|6|, we define the temperature which delineates 
the border between coherent and phase-fluctuating con- 
densates as = L^T/L. Since ni(0), L, and T are 
extracted directly from the images, the definitions of 
and Tfj, relate the coherence properties to experimentally 
measured quantities. The axial momentum distribution 
follows from a Fourier transform of C(s) and is well 
approximated by a Lorentzian of width Ap^ = 
(HWHM), with a — 0.67 for a parabolic no(r) 
The predicted spectral width is therefore aAv^, where 



Av$ = 



2hk L 
2%ML d 



(2) 



In the following we will use a as a free parameter to test 
the theory outlined above. 

Figure shows the measured spectral width Az/jm ver- 
sus A Dispersion of the measured fii(0) and T re- 
sults in the horizontal error bars on Az/^, while verti- 
cal error bars indicate the standard deviation of the fit 
width. The measured widths increase at higher Au^ as 
expected. To compare these data to theory, we need to 
take into account a finite "instrumental" width of the 
Bragg spectra, including the effect of the mirror vibra- 
tions, residual sloshing in the trap, and Fourier broad- 
ening due to the 2-ms pulse length (125 Hz HWHM). 
We assume that all experimental broadenings result in 
a Gaussian apparatus function of half- width wq, to be 
convolved by the Lorentzian momentum profile with a 
half- width aAv$. The convolution, a Voigt profile, has a 
half- width aA^/2 + y/w(, + (aA^) 2 /4. Note that fit- 
ting a Voigt profile instead of a Lorentzian to a spectrum 
gives the same total HWHM to less than 5%, but the 
Lorentzian shape is too predominant to extract reliably 
the Gaussian and the Lorentzian contributions to the 
profile. Using a and wq as free parameters to fit the data 
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FIG. 2: Half-widths at half-maximum Ai/m of the experi- 
mental Bragg spectra versus the parameter Av^ cx h/L^ (see 
text). Vertical error bars are the standard deviations of the fit 
width; the horizontal error bars are the one-sigma statistical 
dispersions of Au^. The solid line is a fit assuming a Voigt 
profile for the spectra. 



of Fig. H we find w G = 176(6) Hz, and a = 0.64(5)(5). 
The first uncertainty quoted for a is the standard devia- 
tion of the fit value. The second results from calibration 
uncertainties on the magnification of the imaging sys- 
tem and on the total atom number, which do not affect 
wq . The agreement of the measured value of a with the 
theoretical value 0.67, to within the 15 % experimental 
uncertainty, confirms quantitatively the temperature de- 
pendence of the momentum width predicted in Ref. |(J. 
The coherence length h/p^ deduced from this measure- 
ment varies between 5.9(8) and 39(4) /im, in the range 
6 < T/T < 36. 

Another important aspect of quasi-condensates, the 
suppression of axial density fluctuations, is investigated 
here through the size of the trapped condensate. In the 
presence of density fluctuations, the resulting interaction 
energy |24j would increase the size with respect to the 
expectation for a smooth density profile. The measured 
axial half-length L after release faithfully reflects the half- 
length in the trap since axial expansion is negligible |25| . 
Figure |21 shows that the measured values of L at various 
temperatures (filled circles) are smaller than the standard 
Thomas-Fermi prediction (diamonds) = 2/j/(M^), 
with fi given by 2fi = hu) (15 N a / 'a) 2 :/ ' 5 |2fJ, where u> = 
(w^) 1 / 3 , a = y/h/(Mu>) and a = 5.32 nm 0. We find 
that this smaller value results from the compression of the 
quasi-condensate by the 3D thermal component (excita- 
tions with energy much larger than Hu>x), and from radial 
quantum pressure. The 3D excited states_contribute neg- 
ligibly to the fluctuations of the phase |6J and only the 
density profile is affected. Using a Hartree-Fock approach 
26, 281, wc find: 



2ff 



n (0) + -j- 
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with the coupling constant g — 4Trh 2 a/M, the ther- 
mal de Broglie wavelength A T = [2nh 2 / (Mk B T)} 1/2 , and 
S3/2( x ) = S^Li x n /n 3 / 2 . The open circles in Fig.|21show 
the solution of Eq. © assuming a parabolic profile, such 
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FIG. 3: Half-length L of the quasi-condensate versus tem- 
perature. Experimental values (• with statistical error bars) 
are compared to three calculations: standard T — Thomas- 
Fermi (♦) and Eq.[3]in the text, with and without radial quan- 
tum pressure (•■£)•• and □ respectively). Calculations use 
the measured temperature and number of condensed atoms, 
as determined by a two-component fit to averaged absorption 
images. Inset: Number of condensed atoms iVo with statis- 
tical dispersion at each temperature. 



that no(0) = 15iVoL~ 3 e 2 /(87r). The aspect ratio e is cal- 
culated according to the theory developed in 29] , which 
takes into account radial quantum pressure. The calcu- 
lated lengths are in agreement with our measurements 
to better than both the statistical error (shown), and 
our estimated calibration uncertainty (not shown) of 4 % 
on the ratio L meSuS J L ca \ c ,. We conclude that a phase- 
fluctuating condensate has the same smooth profile as a 
true condensate, and thus that the axial density fluctu- 
ations are suppressed even when phase fluctuations are 
large. 

In conclusion, we have demonstrated three impor- 
tant features of quasi Bose-Einstein condensates: (i) the 
momentum distribution shape, found Lorentzian; (ii) 
the temperature dependence of the momentum width; 
and (iii) the suppression of density fluctuations. Our 
results are in quantitative agreement with the finite- 
temperature, interactive theory developed in 0, Q sup- 
plemented by a Hartree-Fock treatment of 3D excited 
states. The same method could be applied to investi- 
gate how long range order develops during the conden- 
sate growth (see |l7j and references therein). 
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